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1. Introduction 

The physical behaviour of real gases and liquids is different from that of ideal gases mainly as a consequence of 
the effect of the intermolecular forces. 

Aim of this work is to investigate, in the framework of kinetic theory, the problem of waves propagation in liquids, 
bringing out the role that intermolecular forces play on the behaviour of wave propagation, in particular on the 
dispersion relation. Also, there will be shown the existence of an effect of damping, of the type of Landau damping 
in plasmas, that can be revealed only in the context of microscopic theory, since this effect is strictly related to the 
form of the distribution function and disappears if oscillations are analyzed in the context of macroscopic theory. 
In analyzing waves, one can choose between two methods of approach: either derive the wave equation from 
macroscopic equations (to wit, the fluid-dynamics equations), or in the framework of kinetic theory, i.e., with a 
microscopic approach, as is done in the present work. Here the starting point is the Vlasov equation, particularized 
for liquids; in this equation the effects of interactions between the molecules of the system is accounted for through 
a self-consistent field, that in the present case is presented in Appendix, to which the interested reader is referred. 



* E-mail: domiziano.mostacci@unibo.it 



Wave Propagation And Landau- Type Damping In Liquids 



The Vlasov equation approach is appropriate for liquids, since molecules therein are subjected to simultaneous 
interactions with a large number of surrounding molecules, and hence correlation is negligible. This approach is 
often referred to as mean-field approach (see refs. [1-4]; also [5], pp. 87 and ff.). 

2. Vlasov equation and wave propagation 

The starting point of the present discussion, as mentioned in the introduction, is the Vlasov equation for the 
distribution function /(r, v, t) 

f + (1) 
at or m dv 

In this equation the global effect of molecular interactions is accounted through the Vlasov self-consistent field 
F', which is the average effect in the point considered of the forces from all the surrounding particles, weighted 
on the density distribution of these latter: 

F'(n) = / n(r 2 )^dr 2 (2) 

where f/i i 2(ri,r 2 ) is the pairwise interaction potential between two molecules located in ri and r 2 respectively 
and n(r 2 ) is the number density at r 2 . Given the above definition, the self-consistent field vanishes in any 
homogeneous configuration. 

The Vlasov equation is particularly suitable to study wave propagation in a system where the self-consistent field 
is dominant. In this work, the self-consistent field is derived from a modified Lennard- Jones model (Figure 1), 
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(where r = |ri — r 2 |), as discussed at large in the Appendix, to which the interested reader is referred for further 
details. There, the self-consistent field is calculated for the one-dimensional case that will be considered here, as 

Fl{z) ^ A ^1 (4) 
az 

where the parameter A is given by 



A = 



1 /roN f 
3 V a 
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and is negative for the present case of a liquid, as discussed in the Appendix. Consider now a liquid with no 
external forces (gravity will be neglected here), in some equilibrium described by a distribution function /o(v), 
and a small perturbation defined as follows 
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Figure 1 Modified Lennard- Jones intermolecular potential 



/(«,v,t) = / (v) + ¥ >(*,v,t) (6) 

and hence 

n(z,t) = / f{z,w,t)dav = I / (v)dv+ / tp(z, v,i)dv = n + r)(z,t) (7) 

Possible propagation of the perturbation under these conditions will be, if any, along z the direction. Introducing 
the above prescription into (1) and treating, consistently, the self-consistent field as a perturbation as well, upon 
neglecting terms of order higher than the first, the Vlasov equation becomes 

~xT + Uz ~x — 1 T~ ~x — ~ u y°) 

at oz m az av z 

Taking a Fourier transform from z to fc, and a Laplace transform from t to s, 

~r, \ ., ~,i \ ikANtk, s) dfoifz) - \ 

s<fi(k, v z ,s) + ikv z ip{k, v z , s) -\ = <po(k, v z ) (9) 

m ov z 

where (fio(k,v z ) is the Fourier transform of <p(z,v z ,t) at time t = 0, and ip(k,u z ,s) and N(k,s) are the double 
transforms of the perturbations ip(z, v z ,t) and ri(z,t). 
Rearranging (9) yields 
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ip = 



s + ikv z 

And integration over v z leaves, after rearranging 



<fo(k, v z ) ■ 



ikAN(k, s) df 
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Assuming that the equilibrium distribution function is maxwellian, as is to be expected in the present setting, 
fo(v z ) is calculated as 



foiyz) = J ^ f M (y)dv x dv y = n o(^ 2n x BT 
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and the derivative inside the integral in the denominator of (11) becomes 



df (i 



dv z 



= no\l —2v z /ie 
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where /3 = 2 K b T ■ 

Now the integrals in (11) have to be handled in the complex plane, and the complex extension of the variable v z 
will be noted in the following as us. The above integrals will then be written as path integrals, and (11) rewritten 
as follows 



■i f r cp (k,w)- 



T mJr 



(14) 



m J r dw w — ij 



where the path T is the straight line that lies on the real axis. The function , complex extension of (13), is 

analytic everywhere, and tpo(k, w) can be assumed to be well behaved as well, for physical systems: if this is the 
case, the only singularity in either integrand is the simple pole in w — t£. Without attempting Laplace inverse 
transformation of (14), it is noted that the result in (14) parallels that found by Landau [6], allowance made for 
the factor in front of the integral in the denominator and for the specific form of the equilibrium distribution 
function f (w): the same procedure can be followed, leading to the asymptotic solution [5-7]. At large enough 
times, the evolution is dominated by the rightmost pole of (14), located in the rightmost zero of its denominator, 
call it si, producing asymptotically a simple exponential behaviour: 



fj(k, t) oc exp{sit} 



(15) 



If the real part of Si is negative, it produces a "Landau-type" damping of the propagation. Now, si is a solution 
of the following equation 
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Figure 2 Integration contour for Re[i^] < 

where the result of integration is a complex function of the complex variable s. Equation (16) is the dispersion 
relation of the propagation. To investigate the asymptotic behavior, it is necessary to define an analytic extension 
to the entire complex plane s of the functions resulting from the path integral: this is effected through the choice 
of the contour shown in Figure 2 for the integration path F. 

When Re[s] < , and the pole in w — ij- takes on values located on or below the real axis, the integration path 
is be deformed to circle the singularity, and calculation of the integral yields two terms: the Cauchy principal 
value and the residue in the pole. Hence, it is convenient to look for the solution to (16) separately in the right 
half plane Re[s] > and in the left one Re[s] < 0. The analysis will be conducted with reference to oscillations 
with wavelength such that the wave phase velocity is larger than the mean thermal velocity, that is 
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This hypothesis is introduced on account of the physical consideration that if the mean thermal velocity is larger 
than the phase velocity, the random motion prevails on the orderly motion of the wave and the wave does not 
propagate. 

2.1. Right half plane: Re[s] > 

Calculating the principal value of the integral, (16) becomes 
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where erfc is the complementary error function. In the hypothesis of (17), the er/c function can be conveniently 
expanded as follows [8] 



Wave Propagation And Landau- Type Damping In Liquids 



^ze z \rfc{z) = 1 - -L + A + o 



(19) 



and inserting this expansion into (18) yields 
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By simple inspection, f = is not a root of (20), so the equation can be multiplied throughout by (f ) 4 yielding 
finally 
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The complex variable s may be written as s — oj(a + i): then, in the hypothesis \a\ « 1 
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Introducing these values into (21) 
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Separating real and imaginary parts, the following two equations are obtained 
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From (25) the phase velocity can be calculated: 
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where, again, it should be born in mind that the parameter A is negative, as discussed in the Appendix. Intro- 
ducing it into (26) 
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from which ot = 0, in other words, there is no solution with Re[s] > 0. 
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2.2. Left half plane: Re[s] < 

In this case from (21) and (24) and following the same procedure as in the case Re[s] > 

C(S ' fc) = 1 - ^ + 2^ + — k = (29) 

Introducing the approximation of (22) and (23) and recalling the hypothesis |a| << 1 so that 



eW ^e~^e i2pa ^ (30) 



the real part of the dispersion equation becomes 



(*) + — + 2^ + (fc) ~ ^ (J { 5QC ° S ( 2/3 (fc) "j"""^) "J} (31) 

Now, observing that the expression into parenthesis ~ , the relation between oj and k that gives the waves that 
can propagate in liquids, is the same as (25). For the imaginary part, one obtains 

2An fc 2 /3/c 2 2 \ 4A /33w 5 -2<£ f ( 62u> 2 a\ . ( I32uj 2 a\\ „ 

The expression into parenthesis is « 1 and then the value of a is given by 




(33) 



3. Conclusion 



In the present work, results are obtained from a kinetic theory approach, in particular: the dispersion relation; 
the phase velocity as a function of the intermolecular force; the "Landau- type" damping effect, i.e., a damping 
analogous to that encountered in plasmas, which can also only be seen from a kinetic approach. Starting from 
kinetic theory, some aspects are seen that cancel out when a macroscopic equations approach is followed. 

Appendix 

In the Vlasov equation [7, 9] 

df df F' df 
at or m av 

the effect of molecule interaction is accounted for through a self-consistent field F', to be calculated as 
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F'(r)= / n(r 2 )F 1 , 2 dr 2 (A.2) 
Jv 

where Fi, 2 is the force that a molecule located at position r 2 exerts on the molecule in r and n(r) is the local 
number density at r. The detailed form of the interaction function can be investigated only through quantum 
mechanics and much work has been done in this direction [10, 11]. However the problem is very complex and 
many effects are involved; moreover the structure of the molecules is often not very well known. Therefore the 
existing results contain significant approximations and are applicable only to specific situations. This being the 
case, it becomes essential to resort to a phenomenological potential (7i, 2 . In this work, the following Lennard- 
Jones model, modified to include a distance of closest approach a to account for the non-vanishing dimensions of 
the molecules, will be used to calculate the self-consistent field. The intermolecular potential (7i, 2 (r) (henceforth 
referred to as modified Lennard- Jones model or mLJ) is presented in Figure 1, and is given by 



C/i, 2 (r) = 



ro 



'0 



for r < a < ro 
for r > a 



(A.3) 



To calculate the self-consistent force F' , a molecule located at the point (0, 0, z) will be considered, and the 
force exerted on this by the whole surrounding liquid will be calculated from the mLJ potential. To simplify the 
problem, a system possessing slab symmetry will be assumed, that is one in which density depends only on the 
z - coordinate. Consider then an elementary volume dV at a location defined by the coordinates (r, $,7) in a 
spherical reference system centered in the molecule of interest and with the polar axis along the z direction, see 
Figure 3 [12] 

With the geometry in Figure 3 the force acting on the molecule of interest due to a molecule in (r, ■&, 7) becomes 



Fi. 2 = 4e 



6rg 12rf 



for r > a 



(A.4) 



Now, calling (£, rj, Q the cartesian coordinates of volume dV, the value of r a+1 with a — 12 or a = 6 can be 
calculated as 



ct + l 



(A.5) 



Sines and cosines of the angles in Figure 3 can be expressed in terms of the cartesian coordinates (f , rj, £) 



sin ■& ■- 
cos $ 



v^+^+tc-z) 2 



sin 7 = 



y« 2 +i 2 +(c-z) 2 

The Cartesian component of the force may be rewritten as 



cos 7 ■ 



(A.6) 
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Figure 3 Calculation geometry 
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(A.7) 



(A.8) 



(A.9) 



To obtain the overall force on the reference molecule, integration over the whole volume is performed. It can be 
seen readily that 
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so that there are no x and y components to the force - and this is consistent with the symmetry of the problem. 
As for the z component 



/oo r oo 
- oo J — 



Fi,de*7 = n(0dC(C-*) 



oo J — oo 



oo p oo 



-ooi-oo [£2 +7? 2 + (£_ z )2]- 

Considering the mLJ potential (A. 6), there is a minimum approach distance a 
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F L = 8tt£ 



J-oo (^-C) 5 J z + a(z-0° . 



(A. 12) 



If the density variation is mild, n(z) can be expanded in Taylor series retaining only the first few terms 
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Neglecting terms of order 5 and higher, and substituting into (A. 12), after some algebra the following equation is 
obtained: 



dz 



where the coefficients are given by 
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In the present work, only the first term in (A. 14), i.e., A, will be retained, yielding 



(A. 14) 



(A.15) 
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The question arises of the sign of A, which depends on the value of the ratio: if this is larger than \f?> (i.e., 
approximately 1.20), A becomes negative. Now the value of the distance of closest approach a is dependent on the 
temperature and the density [13], however phenomenological considerations may yield a qualitative answer: liquids 
do not yield easily to compression, and on the other hand have a strong cohesive behaviour; departures from their 
equilibrium density at the given temperature is strongly resisted. This leads one to believe that the coefficient A 
needs to be negative, so that the force opposes the density gradient, trending to restore the equilibrium density. 
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